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CT*. Abstract 

Grover's quantum search algorithm is considered as one of the milestone in the field of 
quantum computing. The algorithm can search for a single match in a database with N 
- records in 0(\/iV) assuming that the item must exist in the database with quadratic speedup 

I over the best known classical algorithm. This review paper discusses the performance of 

'sj" ' Grover's algorithm in case of multiple matches where the problem is expected to be easier. 

Unfortunately, we will find that the algorithm will fail for M > 3N/4, where M is the number 
of matches in the list. 

OO 
O 

1 Introduction 

X. 

. In 1996, Lov Grover [11] presented an algorithm for searching an unstructured hst of N items 
with quadratic speed-up over classical algorithms. His original algorithm targets the case where 
a single match exists within the search space. Much research effort has gone into analysing and 
generalising his algorithm for multiple matches [31 HJ EJ [71 IH] • 

This paper will review the work done by others on solving the unstructured search problem 
on quantum computers as follows: Section 2 provides the general definition of the unstructured 
search problem and some of its applications. Section 3 briefly summarises the work done so far in 
designing algorithms concerning this problem on quantum computers. Section 4 presents Grover's 
algorithm in some detail and the work done by others related to his algorithm, analysing its 
performance and behaviour over the range 1 < M < for both known and unknown number of 
matches M. The paper ends up with a general conclusion in Section 5 about Grover's algorithm. 
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2 Unstructured Search Problem 



Consider an unstructured list L oi N items. For simplicity and without loss of generality we will 
assume that N = 2" for some positive integer n. Suppose the items in the list are labelled with the 
integers {0, 1, A^ — 1}, and consider a function (oracle) / which maps an item i G L to either or 
1 according to some properties this item should satisfy, i.e. f : L ^ {0, 1}. The problem is to find 
any i & L such that f{i) = l assuming that such i exists in the list. In conventional computers, 
solving this problem needs O (N/M) calls to the oracle (query), where M is the number of items 
that satisfy the oracle. 

The unstructured search problem can be considered as a general domain for a wide range of 
applications in computer science, for example: 

• The database searching problem, where we are looking for an item in an unsorted list. 

• The Boolean satisfiability problem, where we have a Boolean expression with n Boolean 
variables and we are looking for any variable assignment that satisfies this expression. 

3 Unstructured Search on Quantum Computers 

Grover's original algorithm exploits quantum parallelism by preparing a uniform superposition 
that represents all the items in the list then iterates both an oracle that marks the desired item by 
applying a phase shift of -1 on that item (e-^\ with 6i = vr) and nothing on the other items (e-^^, 
with 62 = 0) and an operator that performs inversion about the mean (diffusion operator) to am- 
plify the amplitude of the match. The process of this operator includes the operation (2 |0) (0| — I) 
which applies a phase shift of -1 on the states within the superposition (e-"^^, with 0i = vr) except 
the state \0)^^ where it applies nothing (e-'^^, with 02 = 0) (Fig. [2]) [12]. To maintain consistency 
with literature, this operation can also be written as (/ — 2 |0) (0|) which applies a phase shift of 
-1 on the state lO)*^" {e-'^^, with 02 = tt) and nothing on the other states of the superposition (e-"'^^, 
with 01 = 0) together with a global phase shift of -1 (Fig. [3]) [15] . 

It was shown that the required number of iterations is approximately tt/A^/N which is proved 
to be optimal to get the highest probability with the minimum number of iterations [20], if there 
is exactly one match in the search space. 

In [I1[IDI[I21[IS1[IZ], Grover's algorithm is generalised by showing that the uniform superposition 
can be replaced by almost any arbitrary superposition and the phase shifts applied by the oracle 
and the diffusion operator [e-^^ , e-^^ , e-"^^ and e-'^^) can be generalised to deal with the arbitrary 
superposition and/or to increase the probability of success even with a factor increase in the number 
of iterations to still run in 0{^N/M). These give a larger class of algorithms for amplitude 
amplification using variable operators from which Grover's algorithm was shown to be a special 
case. 

In another research direction, work has been done trying to generalise Grover's algorithm with 
a uniform superposition for the case where there are a known number of multiple matches in the 
search space [31 El E], where it was shown that the required number of iterations is approximately 
/A^J N / M for small M/N . The required number of iterations will increase for M > N/2, i.e. the 
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problem will be harder where it might be expected to be easier [12] . Other work has been done for 
a known number of multiple matches with arbitrary superposition and phase shifts [21 HI [TU [161 flS] 
where the same problem for multiple matches occurs. In [HO [18], a hybrid algorithm was presented 
to deal with this problem. It applies Grover's fixed operators algorithm for n / A\J N / M times then 
applies one more iteration using different oracle and diffusion operator by replacing the standard 
phase shifts = vr, = vr with accurately calculated phase shifts 9^ and (l)^ according to the 
knowledge of the number of matches M to get the solution with probability close to certainty. 
Using this algorithm will increase the hardware cost since we have to build one more oracle and 
one more diffusion operator for each particular M . For the sake of practicality, the operators 
should be fixed for any given M and are able to handle the problem with high probability whether 
or not M is known in advance. 

In case of multiple matches, where the number of matches is unknown, an algorithm for esti- 
mating the number of matches (known as quantum counting algorithm) was presented [5l[T8]. In 
[3], another algorithm was presented to find a match even if the number of matches is unknown 
which will be able to work if M lies within the range 1 < M < 3A^/4, otherwise it was suggested 
to use standard sampling techniques. 

Many Grover-like algorithms for solving a wide range of applications have been presented. For 
example, an algorithm for the scheduling problem (intersection problem) [13] and an algorithm for 
minimum finding [9]. 

4 Grover's Quantum Search Algorithm 
4.1 Number of Matches is Known 

In this section, we will present Grover's algorithm for searching a list of size with M matches 
such that 1 < M < A^. We assume that M is known in advance. For our purposes, the analysis 
will concentrate on the behaviour of the algorithm if iterated once, then the behaviour after qc 
iterations. 

4.1.1 Iterating the Algorithm Once 

For a list of size N = 2^, the steps of the algorithm can be understood as follows (its quantum 
circuit is shown in Fig. [T]): 

1- Register Preparation. Prepare a quantum register of n + 1 qubits. The first n qubits all 
in state |0) and the extra qubit in state |1) where it will be used as a workspace for eval- 
uating the oracle Uf. The state of the system Wq'^'^''^ can be written as follows, where 

the subscript number refers to the step within the iteration. [G, 1) in the superscript is the 
diffusion operator used in the algorithm which will be defined later and the iteration number 
respectively: 



^ = |0)«"®|1). (1) 
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Figure 1: Quantum circuit for Grover's algorithm. 

2- Register Initialisation. Apply the Hadamard gate on each of the n + 1 qubits in parallel so 
that the first n qubits will contain the 2"' states representing the list and the extra qubit will 
be in the state (|0) — |1))/ v^, where i is the integer representation of the items in the list: 
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3- Applying the Oracle and Changing Sign. Apply the oracle Uj that gives the amplitudes of 



the matches a phase shift of —1 (e-'^), i.e. Uf \i) 

1 
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Notice that, if /(z) = 0, then |0 © / (z)) = |0) and |l©/(i)) = and if f{i) = 1, then 
|0 © / (z)) = |1) and |1 © / (i)) = |0). Assume that denotes a sum over i which are 
desired matches and Yli " denotes a sum over i which are undesired items in the list. So, the 



system 



(G,l) 



shown in Eqn. [3] can be re- written as follows: 



(G,l) 
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Notice the change of the sign for the states that represent the matches in the search space 
(phase shift of -1), with no change to the state of the extra qubit workspace, which can be 
removed from the system for simplicity. We end-up with a system as follows: 
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Figure 2: Quantum circuit for the diffusion operator G over n qubits. 

4- Inversion about the Mean. Apply the Diffusion Operator G on the first n qubits. The 
diagonal representation of G can take this form (its quantum circuit is as shown in Fig. [2] 
and its quantum circuit with a global phase shift factor of -1 [T9j is as shown in Fig. [3]): 



Q = H®n ^2|o) (0| - J„)i7®", 



(6) 



where the vector |0) used in Eqn. [6] is of length = 2", and /„ is the identity matrix of size 
2" X 2"-. Consider a general system \ip) of n-qubit quantum register: 
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The effect of applying G on produces, 
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j=0 



where, {a) = J2j^=o '^j is the mean of the amplitudes of the states in the superposition, 
i.e. each amplitude a, will be transformed according to the following relation: 
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From Eqn. [5] we can see that there are M states with amplitude — 1/ y/N and N — M states 
with amplitude l/\/iV, so the mean (a) can be calculated as follows: 
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Figure 3: Quantum circuit for the diffusion operator G over n qubits with a global phase shift 
factor of -1 [19]. 



The effect of applying G on the system 
follows: 



Wo^'"^^) shown in Eqn. [5] can be understood as 



a- The M negative sign amplitudes (solutions) will be transformed from —1/^/N to af, 
where af is calculated as follows: Substitute aj = —l/\fN and (a) from Eqn. [TOl in 
Eqn. Owe get: 



_ ( I ^ h _ 2M 

1 (n _ 4A/\ 

W N ) ■ 



b- The (A^ — M) positive sign amplitudes will be transformed from 1/ \/N to 6f , where 
is calculated as follows: Substitute Uj = and (a) from Eqn. [I0]in Eqn. [9] we get: 
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The new system 



1^3*^'^^ ) after applying G can be written as follows, and the mechanism of 



amplifying the amplitudes can be understood as shown in Fig. HI 
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such that, 
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Figure 4: Mechanism of amplitude amplification for Grover's algorithm with N = 8 and M = 1. 

5- Measurement. Measure the first n qubits. The probabilities of the system will be as follows: 

i- Probability Pi^''° to find a match out of the M possible matches can be calculated as 
follows: 



,(l)c 



M{af) 



G\2 



(15) 



9(f)-24(f)Vl6(M)'. 

ii- Probability PnJ^ to find undesired result out of the states can be calculated as follows: 

pWc = (Ar_M)(6ff. (16) 

Notice that, using Eqn. [Ill 

(17) 



^i'^^ + Pi'J"" = M(af )2 + (iV - M){bfY = 1. 



Performance of a Single Iteration 

Considering Eqn. [TTl Eqn. [T^l Eqn. [T3] and Eqn. [T51 we can see that the probability to find 
a solution for a fixed size search space varies according to the number of matches M in the 
superposition (see Fig. [5]). 

From Tab. [1], we can see that the maximum probability of success is always 1.0, and the 
minimum probability (worst case) is always 0.0. The average probability over the possible oracles 
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n, where N = 2^ 


Max. prob. 


Min. prob. 


Avg. prob. 


2 


1.0 


0.0 


0.5 


3 


1.0 


0.0 


0.5 


4 


1.0 


0.0 


0.5 


5 


1.0 


0.0 


0.5 


6 


1.0 


0.0 


0.5 



Table 1: Performance of the first iteration of Grover's algorithm with different size search space. 



for different number of matches is always 0.5. It implies that the average performance of the first 
iteration remains constant even with the increase of the size of the list. 

To verify these results, taking into account that the oracle Uf is taken as a black box, we can 
define the average probability, average{Ps^'*'^), as follows: 



average (pPA = ^ ''CmPs^''^'' 
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where ^Cm = m\{n-m)\ number of possible cases for M matches. We can see that as the 

size of the list increases (A^ oo), average{Ps^^^) shown in Eqn. [T8]remains one-half. 

Classically, we can do a single trial guess to find any match. We may succeed in finding a 
solution with probability pj^'"***'^"') = M/N. The average probability can be calculated as follows: 

N 

average[i^s ) — w ^ '^m^s 

M=l 



_ ^ (19) 

2^ (A/-l)!(7V-Af)! 



1 (N^iy. 



E 

M=l 
1 

2' 



It means that we have an average probability one-half to find or not to find a solution by a 
single random guess, even with the increase in the number of matches, similar to the first iteration 
of Grover's algorithm. 

To compare the performance of the first iteration of Grover's algorithm and the classical guess 
technique. Fig. [5] shows the probability of success of the two algorithms just mentioned as a 
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Figure 5: Probabilities of success 
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^{classical) 



as a function of M/N. 



function of < M/N < 1. 

We can see from Fig. [5] that Grover's algorithm solves the case where M = N/A with certainty. 
The probability of success of Grover's algorithm will be below one-half for M > N/2 and will fail 
with certainty for M = 3N/4. The probability of success of the classical guess technique is always 
over that of Grover's algorithm for M > N/2. 

4.1.2 Iterating the Algorithm 

Before we go further in the analysis of the algorithm after arbitrary number of iterations qc, we 
will re- formulate the equations of the first iteration according to the way used in [3]. 

Initially before the first iteration, we had M states with amplitude = 1/ \/N and N — M 
states with amplitude = 1/\/N . After applying the oracle Uf and the diffusion operator G, the 
new amplitudes af and 6f can be re-written as follows: 
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(20) 
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The iterative version of the algorithm can be summarised as follows: 

1- Prepare a quantum register of n + 1 qubits. The first n qubits all in state |0) and the extra 
qubit in state |1). 

2- Apply the Hadamard gate on each of the n + 1 qubits in parallel. 

3- Iterate the following steps qc times, 

i- Apply the oracle Uf. 

ii- Apply the diffusion operator G on the first n qubits. 

4- Measure the first n qubits to get the result with probability Ps'"^\ 

The system after ^ 1 iterations can be written as follows, 

|w^^^-^> = &?E"i^) + «?E'i^)' (21) 

such that, 

M(aJ^)2 + (iV-M)(6^)2 = l, (22) 

where the amplitudes and 6^ after ^ 1 iterations are defined by the following recurrence 
relations j3]. 



% - % - 

_ N-2M „G , 2(7V-M)lG _ N-2M uG _ 2M ^G 

Solving these recurrence relations, the closed forms can be written as follows [3]: 



(23) 



= sin {{2qG + 1) ^g) , = ^^^W ^^^^^ + ' ^^^^ 



where sin^ {Oq) = M/N and < 6*0' < tt/2. 
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The probabilities of the system will be as follows: 
1- The probability of success after qc iterations is: 



Pi''^) = M {a^Y = sin2((2gG + 1)0g). 
2- The probability of failure after qc iterations is: 



Pif ) = {N-M) {b^f = cos\{2qG + 1)^g). 



(25) 



(26) 



The aim is to find a solution with probability as close as possible to certainty. It was shown in 
[3] that Pns'^^ = when = ~ 26^^) /4:6g, but since the number of iterations must be integer, 
let qc = {tt/4:9g\ where \qG — QgI — -'-/2- And since, sin^ {9g) = M/N, we have for small M/N 
that, Og > sin (Og) = ^M/N, then. 



Qg 
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TT 71 N ^ N 
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(27) 



where [ J is the floor operation. The lower bound of the probability of success using is Pf 
1 — M/N > 0, which is negligible only for small M/N. 



> 




Figure 6: Probability of success of Grover's algorithm using the required number of iterations qc- 

To demonstrate the real behaviour of Grover's algorithm, we may plot the probability of success 
Pg'^^^ using the required number of iterations for any given M. Fig. E] shows this behaviour 
as a function of < M/N < 1. We can see from the plot that the minimum probability that 
Grover's algorithm may reach is approximately 50.0% when M/N = 0.5. The algorithm will 
behave similar to the classical single random guess for M/N > 0.5 since = in that range. 
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For 0.145 < M/N < 0.5, qc = ^ where we can see that for M/N = 0.25 the algorithm will 
succeed with certainty after a single iteration. For M/N < 0.145, > 1 where the algorithm 
will behave more reliably. In an attempt to avoid this drawback in the behaviour for multiple 
matches, it was proposed in [19] that we can double the search space by adding N non-match 
items so that the number of matches will always be less than half the search space and iterate the 



algorithm Tr/A^y2N/M instead of j N jM so it still runs in O y^jN jMy Using this approach 

will increase the space/time requirements to still get the result with probability at least one-half 
when M = N , where we can get the result with certainty in this case if we did not use that 
approach. 

4.2 Number of Matches is Unknown 

In case the number of matches M is unknown, an algorithm that employs Grover's algorithm [3] 
can be used for 1 < M < 3iV/4 which can be summarised as follows: 

1- Start with m = 1 and A = 8/7. (A can take any value between 1 and 4/3) 

2- Pick an integer j between and m — 1 in a uniform random manner. 



3- 



N~l 

Run j iterations of Grover's algorithm on the state: -4^ Yl K)- 



i=0 



4- 



Measure the register and assume i is the output. 



5- 



If /(z) = 1, then we found a solution and exit. 
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Figure 7: The actual behaviour of the functions representing the required number of iterations for 
known (qc = tt/AOg) and unknown (m^ = 1/ sin {26g)) number of solutions, where the number of 
iterations is the flooring of the values (step function). 

It was shown that the total expected number of iterations is approximately Snic ~ A^yN/M 

for small M/N , where ma > 1/ sin (26^0) = O (y'N/M^ for M < 3iV/4. The algorithm works 

only for 1 < M < 3N/A, where for M > 3N/4, a classical sampling techniques can be used. 

The reason that this algorithm will fail for M > 3N/A is that mc is acting as a lower bound for 
qc for M < N/2. It handles the case where go = in a constant manner for N/2 < M < 3N/4. 
However, it will increase exponentially for M > 3N/4 where it is no longer able to approximate qc, 
i.e. using the algorithm in that range means that the expected number of iterations will increase 
exponentially where the problem should be easier, as shown in Fig. [71 

5 Conclusion 

In this paper, we analysed Grover's algorithm over the whole search space. We found that, although 
Grover's algorithm is optimal [20] for a single match in the search space, its reliability may decrease 
for multiple matches, i.e. the behaviour of the algorithm is not reliable over the whole range where 
the minimum probability it may reach is approximately 50.0% when M/N = 0.5. The best 
behaviour is for M/N < 0.145 and in the neighbourhood of M/N = 0.25. The role of Grover's 
algorithm disappears for M/N > 0.5 where the required number of iterations will vanish in that 
range. Grover's algorithm may not be suitable for practical implementation since a practical 
quantum algorithm should be able to handle both the easiest cases and the hardest cases. 
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